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Abstract. We use the Bethe-Salpeter equation in rainbow-ladder truncation to calculate the ground 
state mesons from the chiral limit to bottomonium, with an effective interaction that was previously 
fitted to the chiral condensate and pion decay constant. Our results are in reasonable agreement 
with the data, as are the vector and pseudoscalar decay constants. The meson mass differences tend 
to become constant in the heavy-quark limit. We also present calculations for the pion and rho 
electromagnetic form factors, and for the single-quark form factors of the 77^ and J |'^^■ 
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INTRODUCTION 

Hadrons are color-singlet bound states of quarks, antiquarks, and gluons. Bound states 
appear as poles in the n-point functions of a quantum field theory. Thus a study of the 
poles in the n-point functions of QCD will tell us something about hadrons. 

In the ultraviolet region, these n-point functions can be calculated using perturbation 
theory. For hadronic observables however, we need to understand the nonperturbative, 
infrared behavior of the n-point functions of QCD. The Dyson-Schwinger equations 
[DSEs], which are the equations of motion of a quantum field theory, provide us with a 
tool to study the n-point functions nonperturbatively. For reviews on the DSEs and their 
use in hadron physics, see Li,2«, A,^- 

MESON PHYSICS 



Mesons can be described by solutions of the homogeneous Bethe-Salpeter equation 

f d\ 

^Pout,Pm,P) = J Jj—j^K{pouuPm,kout,km)X{kout,hn;P), (1) 

with Pin, Pout the 4-momenta of the quark and antiquark, subject to momentum con- 
servation: Pin — Pout = r the Bethe-Salpeter amplitude [BSA], and x{koutihn',P) = 
S{kout)T{kouukin',P)S{kin); the kernel K is the qq scattering kernel. This integral equa- 
tion has solutions F at discrete values of = —M^ (in Euclidean metric) of the to- 
tal meson 4-momentum P. Different types of mesons, such as pseudoscalar or vector 
mesons, are characterized by different Dirac structures. The properly normalized BSA 
^{PouuPm,P) completely describes the meson as a qq bound state. 



Since Eq. ([T]) has solutions at discrete values of = —Mf, one does not obtain the 
"complete" spectrum, including the excited states, by solving a matrix equation once; 
instead, one has to repeatedly solve Eq. ([T]) at different values of in order to find the 
mass spectrum. The ground state in any particular spin-flavor channel corresponds to the 
solution with the lowest mass, Mq. Excited states can be found by looking for solutions 
of Eq. © with a larger mass M/ > Mq, and this can indeed be done L^, JJ. 



Rainbow-ladder truncation 

A viable truncation of the infinite set of DSEs has to respect relevant (global) sym- 
metries of QCD such as chiral symmetry, Lorentz invariance, and renormalization group 
invariance. Here we use the so-called rainbow-ladder truncation, in which the qq scat- 
tering kernel is replaced by an effective one-gluon exchange 

where q = po^t — ^out = Pm — hn, and a{q^) is an effective running coupling. The 
corresponding truncation of the quark DSE is 

S{p)-' = ipZ2+mq{^)Z^ + ^j -^Ana{q^)Df^{q)y^S{k)yy, (3) 

where S{p) = Z{p^)/[i / -\-M{p^)] and q = k — p. This truncation is the first term in 
a systematic expansion [8] of the quark-antiquark scattering kernel K; asymptotically, 
it reduces to leading-order perturbation theory. Furthermore, these two truncations are 
mutually consistent in the sense that the combination produces vector and axial- vector 
vertices satisfying their respective Ward identities. 

For the effective interaction we use the 2-parameter model of Ref. [[lo|l 



e2- 1 + (1+^2/^2 ^^)- 



(4) 



with ^{s) = {l-e^')/s, Y,„ = \2/{33-2Nf), and fixed parameters A^^ = 4 and 
Aqcd = 0.234GeV. The remaining parameters, CO = 0.4 GeV and D = 0.93 GeV^, 
were fitted in [[lo|l to reproduce a chiral condensate of (240 MeV)^ and /^r = 131 MeV. 
The first term in Eq. (HI) models the infrared enhancement of the effective qq scattering 
kernel necessary to generate the experimentally observed amount of dynamical chiral 
symmetry breaking [11] . It was introduced in [10] as a finite- width representation of a 
5 -function [jl2|] . which can be interpreted as a regularized 1/p^ singularity in K (v Ml^ . 
The second term ensures the correct perturbative behavior in the ultraviolet region. 



Meson spectroscopy 



In Table [U we give our results for the equal-mass ground states in each spin channel. 
The masses of the light quarks where fitted in [[lo|] to the pion mass (using equal u and 
d quark masses) and to the kaon mass. The light vector and pseudoscalar mesons are 
described very well by this model: not only their masses, but also a wide range of other 
observables agree with experiments, without adjusting any of the parameters, see [4] and 
references therein. Here we apply this model to heavy quarks as well, and use the vector 
mesons J/xjr and T to fix the c and b masses. 

The mass splitting between the pseudoscalar and vector mesons is too large for the 
heavy quarkonium states, but the decay constants are in reasonable agreement with avail- 
able data. On the other hand, the mass splitting between the vector and the scalar mesons 
is too small; and the scalar-pseudoscalar mass difference is reasonable. Also the axi- 
al vector masses are too small, but the mass difference between the scalar and the l"*""*" 
states is in agreement with data, both for the light and for the c and b quarks. Similar 



results for the light quark sector and for the charmonium states were found in Ref. [11711 
with a slightly different model interaction. Presumably corrections beyond ladder trun- 
cation are necessary for the scalar and axialvector masses: there are significant cancel- 
lations between these corrections in the pseudoscalar and vector channels [SJ, but not 
necessarily in the scalar and axialvector channels. 

Over the entire mass range from the chiral limit up to the bottomonium states, the 
pseudoscalar, vector, and scalar masses can be fitted by 

^meson = Q + Ci + C2 , (5) 

where niq is the current quark mass at our renormalization point /i = 19 GeV. The fit 
parameters are Cq = and Ci = 6.94 for the pseudoscalars, Cq = 0.51 and Ci = 7.27 
for the vectors, and Cq = 0.38 and Ci = 8.65 for the scalar mesons, with a common 
parameter C2 ~ 4.6. The fact that the trajectories can all be fitted with (approximately) 
the same value for C2 means that for large masses, the meson mass differences become 
constant: in the limit mq^oa the above fit suggests AM ^ ACi / ^/Oi. Thus this global 
fit indicates that the mass difference My — Mps approaches 0.07 GeV, whereas Ms — Mps 
approaches 0.4 GeV for heavy quarks; our numerical results however do not exclude that 



TABLE 1. Masses and leptonic decay constants for equal-mass ground state mesons. Exper- 
imental data are from Ref. 1 15], with the exception of f-q^ lfl6ll . 



quark flavor 


Mps 


/ps 


Mv 


/v 


M(0++) 


M(l+-) 


M(l++) 


up/down 


0.1385 


0.131 


0.743 


0.207 


0.672 


0.83 


0.91 


expt. 


0.135,0.140 


0.131 


0.775 


0.221 


0.985 


1.23 


1.23 


strange 


0.697 


0.183 


1.076 


0.260 


1.081 


1.17 


1.25 


expt. 






1.020 


0.229 








charm 


2.908 


0.381 


3.098 


0.421 


3.250 


3.26 


3.33 


expt. 


2.980 


0.335 ±0.075 


3.097 


0.416 


3.415 




3.51 


bottom 


9.38 


0.66 


9.46 


0.62 


9.72 


9.73 


9.75 


expt. 


9.30 




9.46 


0.715 


9.86 




9.89 




FIGURE 1. Meson masses (left) and mass differences (right) as function of current quark mass, 
normalized to the up and down quark masses. The vertical dashed lines indicate physical quark masses. 



the coefficients Ci are identical for the pseudoscalar and vector mesons, in which case 
this mass difference vanishes in the heavy quark limit. 

This is indeed consistent if we look at the actual mass differences we find, see the 
right panel of Fig. [T] the mass difference My — Mps decreases with increasing quark 
mass, it is about /S.m ~ 0.06 GeV for at 2m^, and still decreasing. Similarly, the mass 
difference May — ^^s appears to vanish in the heavy quark limit, but the differences 
Ms — Mps and May — M^s clearly remain nonzero and appear to go to a constant of 
about AM 0.35 GeV. However, one should keep in mind that the model was fitted to 
the pion decay constant and the chiral condensate; implicitly we may have incorporated 
corrections beyond the ladder kernel in our model for the effective qq scattering kernel. 
Higher-order corrections affect light quarks differently than heavy quarks 

The corresponding quark mass functions are shown in Fig. [21 and summarized in 
Tabled Our current quark masses are in good agreement with conventional values [ Tsll 
of both the light and the heavy quark masses. For the light quarks, the nonperturbative 
mass function Mq{p^) is significantly larger than the perturbative quark mass mq{ii) 
at p = 2 = /i, indicating that chiral symmetry breaking sets in well above this scale. 
The momentum dependence of Mc^b{p^) is much less dramatic. Nevertheless, there is a 



TABLE 2. Current quark masses mq{pL) at /i = 19 GeV, scaled down lo pL — 2 GeV and to ji = rriq 
using one-loop pQCD, together with the dynamical mass function M{p^) at several values of p^. 



m^(19) mq{2) 




M,(p2=M,(/72)2) 








chiral limit 




0.392 


0.010 


0.477 


0.594 


0.0037 0.005 




0.401 


0.017 


0.499 


0.610 


0.0838 0.118 




0.556 


0.168 


0.689 


0.845 


0.827 1.17 


1.30 


1.42 


1.31 


1.61 


2.00 


3.68 5.65 


4.46 


4.30 


4.46 


4.52 


5.33 
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FIGURE 2. Dynamical quark mass function Mq{p^) for u = d, s, c, b and chiral quarks. 



significant difference between the dynamical mass in the region relevant for qq bound 
states, namely ~ ~5^meson the timelike region, and Mc^biP^) in the spacelike 
region, even for b quarks. For 0< < — \M^^^g^, the mass function of the heavy quarks 
is in fact quite close to the typical pole masses used in non-relativistic calculations of 
charmonium, mc°^^ ~ 1.47 to 1.83 GeV and bottomonium, in^°^^ ^ 4.7 to 5.0 GeV [|15|] . 

Electromagnetic form factors 

The qqy vertex is the solution of the renormalized inhomogeneous Bethe-Salpeter 
equation with the same kernel K as Eq. ©. Thus for photon momentum Q, we have 

f d'^k 

^^{PouuPin) =Z2Y^l+ J j2^KiPouuPm,Kuuhn) S{kout)T^{kouukm)S{kin) , (6) 

with 7?out and the outgoing and incoming quark momenta, respectively, and similarly 
for ^out and ^in, with pout — Pin = ^out — ^in = Q- The ladder truncation for Eq. Q, 
in combination with the rainbow truncation for the quark propagators and impulse 
approximation for electromagnetic form factors, satisfies the vector Ward-Takahashi 
identity and electromagnetic current conservation is guaranteed. 

Also note that solutions of the homogeneous version of Eq. ^ define vector meson 
bound states with masses My = —Q^ at discrete timelike momenta Q^. It follows that 
has poles at those locations. Thus the effects of intermediate vector meson states 
on electromagnetic processes can be unambiguously incorporated by using the properly 
dressed qqy vertex rather than the bare vertex 7^ ii sll . 

Consider for example the 3-point function describing the coupling of a photon with 
momentum Q to the quark a of a meson ah, with initial and final momenta P±jQ 

Al{P,Q) = iN,j^Tr[Tl{q^,q^)x''\q^.q)S'{qr'x^''{q.q-)], (7) 




strange quark 

r up and down quarks 

chiral limit 



TABLE 3. Static electromagnetic properties of pseudoscalar and 
vector mesons {n and p) and cc mesons (fictitious). 







'V.E 




'V,M 




up/down 


0.44 


0.54 


2.01 


0.49 


-0.41 


charm 
lattice [23] 


0.048 
0.063(1) 


0.052 
0.066(2) 


2.13 
2.10(3) 


0.047 


-0.28 
-0.23(2) 




FIGURE 3. Single-quark form factors: n and p (left) and cc pseudoscalar and vector mesons (right). 



with q = k — jP and q± = k + jP± jQ. The corresponding single-quark elastic form 
factor F" of a pseudoscalar meson is defined by 

2P^F%Q^) = Al{P,Q). (8) 

Vector mesons have three elastic form factors, commonly referred to as the electric, 
magnetic, and quadrupole form factors Ge(2^), Gm(2^), and Gq{Q^). The electric 
monopole moment (i.e. the electric charge), magnetic dipole moment and the electric 
quadrupole moment follow from the values of these form factors in the limit ^ — > 0: 
Ge(0) = 1 (constrained by current conservation), Gm(0) = /i, and Gq(0) = .3. 



Our results for the pion form factor [|18L 11911 are in good agreement with the data, 



both in the spacelike region [20] and in the timelike region; the charge radius agrees 
very well with the experimental value (r|) = 0.44 ±0.01 fm^ [21], see Table [31 The 
vector charge radius 1221] is slightly larger than the pseudoscalar radius, both for light 
quarks and for charm quarks. This suggests that the vector states are broader than the 
corresponding pseudoscalar states, assuming that the charge distribution is indicative of 
the physical size of the bound state. This agrees with the naive intuition that a more 
tightly bound state is more compact than a heavier state with the same constituents. 
For charm quarks this difference is significantly smaller than for up and down quarks, 
in agreement with recent lattice calculations [i23i] . The magnetic moment appears to 
be surprisingly independent of the quark mass; the quadrupole moment decreases with 



increasing quark mass ^22^. Recent lattice simulations (l^ agree quite well with our 
results for the moments of the single-quark form factors of the 7/ 1//. 

In Fig. [3] we see that both the pseudoscalar Fj^ and all three vector form factors 
diverge in the timelike region as —0.55 GeV^, corresponding to the vector-meson 

poles in the dressed quark photon vertex. Similarly, the single-quark form factors of the 
rjc and 7/i// diverge as 2^ ^ —9.5 GeV^. However, it is only the pion form factor that 
can be described by a vector meson dominance [VMD] curve, ^ Mp / [Q^ +Mp], over 

the entire 2^-region shown. The p form factors drop significantly faster [22] than a 
VMD curve, as do the cc form factors, both for pseudoscalar and vector states. 
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